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This Letter theoretically discusses quasiparticle states and nuclear spin relaxation rates T^ 1 
in a quasi-one-dimensional superconductor (TMTSF^PFe under a magnetic field applied parallel 
to the conduction chains. We study the effects of Josephson-type vortices on by solving the 
Bogoliubov de Gennes equation for p-, d- or /-wave pairing interactions. In the presence of line nodes 
in pairing functions, Tj - is proportional to T in sufficiently low temperatures because quasiparticles 
induced by vortices at the Fermi energy relax spins. We also try to identify the pairing symmetry 
of (TMTSF) 2 PF 6 . 

PACS numbers: 74.70.Kn, 74.20.Rp, 74.25.Op, 76.60.-k 



Recently, much attention has been focused on su- 
perconductivity in a Bechgaard salt (TMTSF) 2 PF 6 1], 
where a superconducting phase appears next to a spin 
density wave phase in the pressure-temperature phase 
diagram Q. So far a number of theories have pro- 
posed unconventional superconductivity in this quasi- 
one-dimensional (Q1D) organic superconductor: spin- 
singlet ci- wave d, 0, H 0|, triplet p-wave 
triplet /-wave symmetries [T3, EH EH EH 
though pairing interaction along the chain direction is 
a common conclusion among these theories, the pair- 
ing symmetry itself has been controversial. Experimen- 
tally, the spin-triplet pairing has been suggested from 
the unchanged Knight shift across the superconducting 
transition temperature T c in the NMR experiment [161 ] 
and the large enhancement of H C 2 exceeding the Pauli 
limit [171 ]. The nuclear spin relaxation experiment is a 
powerful tool to identify the orbital part of the pairing 
function. The relaxation rate T-j - should exhibit an ex- 
ponential temperature(T) dependence for fully gapped 
superconductors, while it is proportional to T 3 in the 
presence of line nodes in the gap. Such a theoretical anal- 
ysis was applied to the TMTSF salts by Hasegawa and 
Fukuyama 0Jn (TMTSF) 2 C10 4 , the relation T 1 ~ 1 ocT 3 
as well as the absence of the coherence peak were ob- 
served at zero ma gne tic field, which indicates the pres- 
ence of line nodes [Hj]. In (TMTSF) 2 PF 6 , however, Lee 
et al. [l|| observed T-j -1 oc T in low temperatures under 
a magnetic field at 1.43 T. Under a magnetic field, quasi- 
particles induced at the Fermi energy may cause the spin 
relaxation. In fact, two of the present authors and other 
groups have shown that the relation T~ 1 cx.T holds in a 
number of nodal superconductors [20L HH [HI, HH • In the 
NMR experiment by Lee et al. [l6|, however, the mag- 
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FIG. 1: (a) Fermi surface in the Q1D tight-binding model 
(£7(k)=0) in the k a -kb plane. k c dependence of the Fermi 
surface is negligible in this scale of the figure, (b) The gap 
function Aoip(k a ) along the Fermi surface for p-, d- and /- 
wave pairing symmetries at T=0 and H—Q. 



netic field is applied parallel to the conduction chains. At 
first sight, this experimental configuration seems to re- 
move such quasiparticle contributions from T-j - because 
the order parameter is mainly localized around the chain, 
and the vortices pass through the open space between 
the chains. Thus order parameters along the chain are 
expected to be unchanged from their value at the zero 
magnetic field 0, HE HH • However, no theoretical study 
has ever tried to reveal effects of such Josephson vortices 
on quasiparticle states and Tj - . 

There are two main purposes to this Letter: (i) study- 
ing temperature dependences of Tj~ in the Q1D organic 
superconductor in the presence of vortices, and (ii) identi- 
fying the pairing symmetry of (TMTSF^PFg. For these 
purposes, we calculate quasiparticle states and order pa- 
rameters self-consistently in the presence of the magnetic 
field along the chain direction based on the Bogoliubov- 
de Gennes (BdG) equation. Three pairing symmetries 
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(i.e., p-, d- and /-wave) are assumed on a tight-binding 
lattice. We will show that the relation T 1 _1 ocT holds at 
low temperatures for d- and /-wave symmetries. 

We consider a three-dimensional tight-binding model, 
where the transfer integrals between the nearest neigh- 
bor sites are — t a , — 1& and — t c in the a, b and c di- 
rections, respectively. The band dispersion is given by 
E(k)—— 2t a cosk a — 2th cos kb— 2t c cosk c — \i. To describe 
the electronic structure of (TMTSF) 2 PF 6 , we set t a : 
tf, : t c =10 : 1 : 0.03 and tune the chemical potential /i 
to keep the band filling at quarter filling (i.e., (n)=0.5). 
The Q1D Fermi surface for these parameters is shown in 
Fig- [Ha). We introduce a spin-singlet pairing interaction 
between the second nearest neighbor sites in the chain di- 
rection (a axis) in the case of the d-wave symmetry [13] ■ 
For the p-wave (/-wave) pairing symmetry, we introduce 
a spin-triplet pairing interaction [28[ between the second 
(fourth) nearest sites in the a direction. The pairing func- 
tions result in ip(k a )— cos 2k a , sin2fc a and sin4fc Q for the 
d-, p- and /-wave pairings, respectively. By applying the 
Fourier transformation in the a axis, the Hamiltonian is 
given by 

"H = Kk a ,i,j a j,<r a i,e 
k a ,j',i,(7 

+ ^{AfV#aH,iai,T + ^MkaH^alj, (1) 

where, Kk a ,ij= - U,j + 5 it j(-2t a cosk a - fJt), a\ t(T {ai,a) 
is a creation (annihilation) operator of an electron with 
spin a(— t or j), an d i—(ib,ic) is a site index in the be 
plane. For a magnetic field applied along the a direc- 
tion \H.=(H, 0, 0)], the transfer integral in the be plane 
becomes ?jj=tjjexp[i^ f* J A(r) • dr], where </>o is the 
flux quantum and the vector potential is taken to be 
A(r) = |Hxr in the symmetric gauge. 

From Eq. fTJ) , the BdG equation is derived for each k a 

Ms )■<*> 

where u e> i and v e ^ are the wave functions belonging to 
the eigenvalue E e , and Dk a ,ij=ip(k a )Ai5ij. The self- 
consistent equation for the pair potentials is given by 

A^U^u^vl^k^fiEe), (3) 

e 

where f(E) is the Fermi distribution function and U is 
the strength of the pairing interaction. 

We consider that two vortices accommodate in a unit 
cell of 20x6 lattice sites in the be plane and form a tri- 
angular lattice. Note that the coherence length becomes 
anisotropic in the be plane, (i.e., : £ c oci^ 2 : t^ 2 ). By 
introducing the quasimomentum of the magnetic Bloch 
state, the wave functions are calculated for sufficient 
quantities of unit cells connected by the periodic bound- 
ary condition [2(|. We iterate the calculation of Eqs. @ 



and ([3]) alternately until a self-consistent solution is ob- 
tained. Using the self-consistent solution (u e ,i, Ve,i an d 
E e ), the local density of states (LDOS) is given by 

N(E, Ti )= J2{K,i\ 2 t( E " &e) + WA 2 5{E + E e )}. (4) 

The screening current J(r) is also calculated. For in- 
stance, the current component in the b direction is given 
by J b (r0odm[i M+g £J U£ , lU * i+g /(£; £ ) + u e > M+& (l - 

f(E e ))}], where i + b denotes the nearest neighbor site 
to i in the 6-direction [2(|. In the calculation, we take 
U=— 38tb (giving Ao=2.4t{, and T c =1.05tb) for d-wave 
pairing, U——15tb (Ao=0.62t& and T c =0.35tb) forp-wave 
pairing, and U=—27.5tb (Ao—1.3tb and T c —0.55tb) for 
/-wave pairing with Ao=Aj(iJ=0). Here U is changed 
depending on pairing symmetries so that we have similar 
gap values in the three symmetries. The ratio Ao/T c de- 
pends on the anisotropy of |^i(fc )|. FigureQJb) shows the 
superconducting gap Aoip(k a ) along the Fermi surface at 
H =0 for each symmetry. 

First we briefly summarize the density of states (DOS) 
spectra N(E, r,-) in the absence of the magnetic field (i.e., 
in the uniform state) as shown by the dashed lines in 
Fig.d The dotted lines in Fig. [2] denote the DOS in the 
normal state. In the p-wave symmetry, no nodes exist in 
the gap along the Fermi surface as shown in Fig. QXb). 
Correspondingly, the full gap structure can be seen in 
the DOS as shown in Fig.^a). Two peaks at E—±0.6Mb 
reflect the large enhancement of the DOS at the gap edge. 
In the d- and /-wave symmetries, the gap functions have 
line nodes at /c a =±7r/4 on the Fermi surface. As a result, 
DOS shows a V-shaped gap structure at H=0, (i.e., E- 
linear DOS at low energies) as shown by the dashed lines 
in Figs. (2Kb) and 12(c). 

Before turning to the DOS in the presence of vortices, 
a typical spatial profile of the pair amplitude under the 
magnetic field should be discussed. In Fig. [3(a), we show 
| Ai| for the d-wave symmetry for instance. Since the 
vortex center is located in the space between conduc- 
tion chains, |Aj| is suppressed very slightly even near 
the vortex. The screening current flows around the vor- 
tex as shown in Fig. [3(b), where the 6-component of 
the supercurrent is presented. Compared to the con- 
ventional vortex structure, the screening current flows 
in much wider region around the vortex. These features 
are typical characters of the Josephson vortex in layered 
superconductors 0, [H, HI] . 

Second, we show the LDOS in the vortex state as 
presented by the solid lines in Fig. [21 where we show 
N(E, r) at the lattice site nearest to the vortex. In 
the fully-gapped p-wave symmetry, quasiparticles have 
bound states around the Josephson-type vortex with the 
energy near the full gap energy [25| . The two smaller 
peaks in LDOS around E^±0.56tb reflect such bound 
states. The larger peaks around E^±0.73tb correspond 
to gap edges smeared by the magnetic field. In the p- 
wave, there is no quasiparticle states at the Fermi energy 
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FIG. 2: The LDOS N(E, r) at the site nearest to the vortex core is presented at T—0 for the p-wave in (a), d-wave in (b) and 
/-wave in (c). The LDOS in the vortex state is plotted with solid lines. The results in the uniform state at zero field and those 
in the normal state are shown with the dashed and dotted lines, respectively. 




FIG. 3: Spatial profile of the pair amplitudes and supercur- 
rents in a unit cell are shown for the d-wave symmetry at 
T—0. In (a), the amplitudes of the pair potentials (|Ai|) are 
presented. In (b), the amplitudes of the b component of the 
screening current (J(,(r)) are presented. 

and the wide gap remains even in the presence of vor- 
tices. Such behaviors are typically found in the quasi- 
particle states around the Josephson vortex in full gap 
superconductors. In the d-wave [Fig. H{b)] and the f- 
wave [Fig. [2c)] symmetries, on the other hand, V-shaped 
superconducting gap is filled with low energy quasipar- 
ticles induced by vortices. Although suppression of A, 
around the Josephson vortex is very small as shown in 
Fig. [3ja), vortices surely induce quasiparticles at the 
Fermi energy. Quasiparticles cannot form a zero-energy 
peak of the vortex core state because of the small sup- 
pression of Aj [20l.l2ll]. In both d- and /-wave, there is no 
bound states of quasiparticle in vortices because of line 
nodes. The calculated results in Fig. [5] show characteris- 
tic low energy quasiparticle excitations depending on the 
gap functions under the magnetic field. 

Next we discuss the T-dependence of T^ 1 on the basis 
of the calculated results in Fig. [2 The spin-spin correla- 
tion function x±( r i> v v > i^n) can be calculated from the 
wave functions of the BdG equation. The nuclear spin 
relaxation rate is given by 0, IMl 

RfaTi,) = Imx+,_(ri,r^,in n ^fi + i77)/(O/T)| n ^ 

X7rTf'(E e )5(E e -E e ,). (5) 
We choose rj=iv because the nuclear spin relaxation at a 



local lattice site is dominant. A relation 5(x)=Tr~ 1 lm(x — 
irf) is used to handle the discrete energy levels of the finite 
size calculation. We typically use 77=0.02^. In Fig. [4] we 
show T^ 1 as a function of T at the lattice site nearest 
to the vortex core with solid lines. Although the relax- 
ation time defined by Ti(r) = l/i?(r, r) depends on r, re- 
sults in the following are almost independent of r. In 
the normal state, we reproduce T-linear behavior of the 
Korringa law as shown by the dotted lines in Fig. [U In 
the p-wave [Fig. H[a)], T-f 1 shows an exponential T de- 
pendence for both the vortex and uniform states. Effects 
of the vortices on 1 is almost negligible in low temper- 
atures because there are no quasiparticle states around 
the Fermi energy even in the presence of magnetic field as 
shown in Fig.[2t[a). By contrast, in the <i-wave symmetry 
in Fig. [DJb) and the /-wave symmetry in Fig. 2{c), T^ 1 
deviates from T 3 -behavior and becomes proportional to 
T at low temperatures as observed in the experiment flf3| 
(see the double logarithmic chart in the insets). This 
implies that the quasiparticles at the zero-energy relax 
spins as they do in the normal state. Although the tem- 
perature at which crosses over from cx T 3 to T seems 
to be somewhat below T c compared to the experimental 
result, [l6| this is mainly because we have taken Ao / i& and 
T c /tf, to be larger than in the actual material due to the 
restriction in the numerical calculations. We have con- 
firmed that this crossover temperature becomes higher 
for smaller Ao, so that a more quantitative agreement 
with the experiment is expected for more realistic values 
of A /t b . 

Finally we discuss the pairing symmetry in the actual 
TMTSF superconductor. The ma gnet ic field is applied 
along the a axis in the experiment [l6| to remove the ef- 
fect of the vortices on T^. However our results show 
that the effect of the vortices is not negligible in un- 
conventional superconductors with line nodes and that 
T^ 1 is proportional to T at sufficiently low temperatures. 
When large magnetic field is applied in other directions, 
the effect of vortices can be even more remarkable. In 
fact, T 1 ~ 1 ocT was also observed under a magnetic field 
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FIG. 4: The nuclear spin relaxation rate T^ 1 is shown as a function of temperature for (a) p-wave , (b) d-wave, and (c) /-wave 
symmetries. Solid lines denote T^ 1 at the lattice site nearest to the vortex core. The results for the uniform state at H—0 
and those in the normal state are plotted with dashed and dotted lines, respectively. Insets in (b) and (c) show the double 
logarithmic chart of the results. 



in the b direction [29| . In a very weak magnetic field, 
T^ 1 is considered to be proportional to T 3 as observed 
in the experiment for (TMTSF) 2 C10 4 This further 
indicates that (i) the gap has line nodes that intersect 
the Fermi surface and (ii) the T-linear behavior at high 
magnetic field is caused by the effect of the vortices. We 
cannot discriminate between d- and /-waves solely from 
Tf , but if we combine this with the Knight shift [16| and 
the H C 2 measurement [l7j . triplet /-wave pairing seems 
to be the most promising candidate for the pairing state. 

In summary, we have studied the quasiparticle states 
and the temperature dependence of the nuclear spin re- 
laxation rate T^ 1 in a quasi-one-dimensional organic su- 
perconductor (TMTSF) 2 PF 6 . We have considered the 
situation in which a magnetic field is applied along the 
chain direction as in the experiment. As for the super- 
conducting state, we have assumed three pairing symme- 
tries (p-, d- and /-wave). Although the suppression of 
the order parameters near the Josephson vortex is very 
small, the Josephson vortex induces quasiparticles at the 
Fermi energy when the gap has line nodes (d- and /- 



wave symmetries). As a result, the LDOS at the Fermi 
energy near the vortex becomes finite in the magnetic 
field. Since T^ 1 at low temperatures is very sensitive to 
the LDOS near the Fermi energy, T^ 1 is proportional to 
T in d- and /-wave symmetries. On the other hand in the 
fully gapped case (p-wave symmetry), no quasiparticles 
is induced at the Fermi energy even in the magnetic field, 
so that the exponential dependence of T^ 1 on T remains 
unchanged. We conclude that the relation T^ 1 oc T ob- 
served in (TMTSF)2PF6 [16| is an evidence for line nodes 
in the gap. Together with several experimental observa- 
tion indicating spin-triplet pairing [ill [Tt| . our results 
suggest that the /-wave symmetry is the most promising 
candidate for Q1D organic superconductors. 
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